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From the discovery of the neutrino to the precision neutrino oscillation measurements in KamLAND, nuclear 
reactors have proven to be an important source of antineutrinos. As their power and our knowledge of neutrino 
physics has increased, more sensitive measurements have become possible. The next generation of reactor 
antineutrino experiments require more detailed simulations of the reactor core. Many of the reactor simulation 
codes are proprietary which makes detailed studies difficult. Here we present the results of the open source 
DRAGON code and compare it to other industry standards for reactor modeling. We use published data from 
the Takahama reactor to determine the quality of the simulations. The propagation of the uncertainty to the 
antineutrino flux is also discussed. 
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1. Introduction 

Reactors are a copious source of antineutrinos, and there is long history of exploiting them for measurements 
of fundamental neutrino properties. The first reactor antineutrino experiment w as in fact the discov ery of 



or lundamentai neutrino prop i 
the neutrino by iReines et al.l 



1960j |. More recently, the KamLAND expe riment, Gando et al.l 2011 1, used 



the reactors of Japan to prove that the effects seen in solar experiments, lAharmim et al. 2011 1, Abe et ahl 
(201lj . lBelliniet al.l (201 ij . is due to a non-zero neutrino mass and the resulting oscillation phenomenon. To 
complete our understanding of neutrino oscillation, several experiments are using reactors to measure the only 
unknown mixing ang le ^13. These experiments include Double C hooz, the experi ment that motiva tes this work, 
lArdellier et al.l|2006| . as well as the Daya Bav. lGuo et al.l [2007| . and the RENO lAhn et al.l I2OIOI exp eriments. 
Reactors are also good sou rces for experi ments searching for sterile neutrinos, iBowden e t al. coherent 
nuclear scattering, Hagma nn and Bernstei n 2004J, and the neutrino magnetic moment, lOarakchcva 2004] . In 
addition, there is al so a large effort to use antineutrinos to monitor reactors for nonproliferation purposes, 
iBowden et al.l j2007| . All o f these efforts r equire simulations for their particular reactor cores, and therefore 
need a code like DRAGON, iMarleau et~aLl [l994]. 

Antineutrinos are produced when fission fragments /3-decay in the reactor core. There are four primary fissile 
isotopes that produce antineutrinos in the energy range most relevant to experiment. They are ^'^^U, ^'^^U, 
^■^^Pu, and ^''^Pu, and they produce neutrinos with energies < 8.5 MeV. The flux of antineutrinos that arrives 
at a detector can be expressed as 



j9 isotopes 
d 71,-, 



dE,-At 



f^{t)S^{E,) 



(1) 



where fi{t) is the fission rate of one of the four primary fissile isotopes and Si{Ey) is the spectrum of neutrinos 
_ ■ emitted per fission of that isotope. These spectra a r e extracted from /3- s pectr u m measurements a t the ILL 
k> ', research reactor bv ISchreckenbach et al. [l981 , 1985 1. Ivon Feilitzsch et al.l [1982| |. iHahn et all 1989| |. It is the 
^ ■ re-analysis of these s pectra by Mueller et al.l 201 1| and Huber 201 1| that has led to the "Reactor Anomaly" , 
\ iMention et al.l 2011 1. The DRAGON simulation provides the fission rates, fi{t), based on a detailed description 
of the reactor design and operating conditions during a given time period. 

To understand the accuracy of the DRAGON predict ions, it is necessary to compare the simulation results 
to data from nuclear reactors. In iNakahara et ahl [2001 1. the Japanese Atomic Energy Research Institute has 
published the mass inventories of the key isotopes from spent fuel rods used in the Takahama-3 reactor core. 
This data has been used to benchmark many widely-used simulation codes. We can use this data to compare 
dragon's results to those from these other codes, a nd to understand t he systematic uncertainties in these 
results. This work is a subset of the effort presented bv I Jones et al. 1 12011| . 



2. Reactor Basics 

There are two main types of nuclear reactors: boiling water reactors (BWR) and pressurized water reactors 
(PWR). Most modern reactors, including Takahama-3, are PWRs. The cores of such reactors are formed by 
fuel assemblies. Assemblies in turn are made up of fuel rods arranged in a grid. In this grid, some locations 
are used for instrumentation and some for shaping the neutron flux with special gadolinium rods. A fuel rod is 
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Table I: The position of samples taken from rod SF97, with the corresponding moderator temperatures, and burnup 
values for that sample. Measurements are in mm from the top of the rod. The bottom of the rod is at 3863 mm. 



Sample 


Position 


Moderator Temp. 


Burnup 




[mm] 


[K] 


[GW-days/ton] 


1 


163 


593.1 


17.69 


2 


350 


592.8 


30.73 


3 


627 


591.5 


42.16 


4 


1839 


575.8 


47.03 


5 


2926 


559.1 


47.25 


6 


3556 


554.2 


40.79 



typically 4 m long, and fashioned out of tiny cylindrical fuel pellets measuring 1 cm in diameter and 1 cm in 
height. The fuel pellet is composed of UO2 enriched to a few percent in ^ssy |-,y weight. The structure of the 
rod is maintained by Zircaloy cladding. Zircaloy is a zirconium alloy chosen for its high melting point. 

The assembly is the elementary unit of fuel in the core. It arrives at the power plant assembled with fresh 
fuel, ready to be installed. Reactors refuel approximately once a year. The time between refueling is called a 
fuel cycle, and during each fuel cycle approximately a third of the fuel assemblies will be fresh. During each 
refueling, the assemblies are arranged very precisely to create an approximately uniform neutron flux across the 
core while extracting the most power out of the fuel. 



3. DRAGON 

Since the fuel assembly is the building block of the reactor core, it is logical that DRAGON simulates individual 
assemblies. To model the whole reactor core with DRAGON, one needs to assume that the spill- in of neutrons 
from neighboring assemblies matches the loss of neutrons from spill-out. The results presented here indicate 
that this is good assumption. For those interested in modeli ng the full c o re mo re accurately, DRAGON results 
can be interfaced with the full core simulation DONJON, IVarin et al. I (2004j . The DRAGON code is open 



source^ making it attractive for use in antineutrino studies. In fact, small modifications are made in DRAGON 
for this work to allow the extraction of the fission rates for the antineutrino flux calculation. 

Codes that simulate reactors must solve the neutron transport equations. The deterministic approach uses 
simplifications to allow the direct solution of these equations. In contrast, Monte Carlo based codes use the 
random generation of a large neutron sample to solve the equations statistically. DRAGON is a deterministic 
2D lattice code. This means that it simulates the assembly by solving the neutron transport equation explicitly 
on a 2D lattice of fuel cells roughly corresponding to the fuel rods. For a given time step in the evolution of 
an assembly, the neutron transport equation is solved, and the solution is used to evolve the fuel composition 
according to the Bateman equations. The deterministic technique is fast compared to Monte Carlo based 
approaches. Symmetries in the assembly geometry can be exploited to further simplify and accelerate these 
calculations. 



4. The Takahama Benchmark 

Takaham a-3 is PWR reactor loc ated in Japan. Three fuel rods in two fresh assemblies were part of the study 
outlined in iNakahara et ahl 2001 1. We focus on fuel rod SF97 because it has the largest exposure or burnup, 



and any cumulative systematic effects will be maximized. The assembly containing SF97 was present in three 
consecutive fuel cycles of 385, 402, and 406 days with 88 days and 62 days of cool-down time between cycles. At 
the end of the exposure, six sample disks, 0.5 mm in width, were removed at the positions indicated in Table HI 
The samples were dissolved and a chemical separation was performed. Mass spectroscopy was used to determine 



^http: / /www.polymtl.ca/nucIeaire/DRAGON/en/index.php 
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Figure 1: The burnup of the six SF97 samples as determined using the Nd technique. 



Table II: Key parameters for the DRAGON simulation. 



Parameter 


Value 


Moderator Density 


0.72 g/cm^ 


Moderator Temperature 


600.0 K 


Cladding Temperature 


600.0 K 


Fuel Temperature 


900.0 K 


Fuel Density 


10.07 g/cm^ 


Fuel Cell Mesh 


1.265 cm 


Fuel Rod Radius 


0.4025 cm 


Fuel Cladding Radius 


0.475 cm 


Guide Tube Inner Radius 


0.573 cm 


Guide Tube Outer Radius 


0.613 cm 



the mass inventories of these samples. For the four isotopes of greatest interest to us, 235^^ ^^8^^ ^-^^Pu, and 
^"^^Pu, the uncertainty is <0.1% for uranium isotopes and <0.3% for plutonium isotopes. iNakahara et al.l 200l| . 

Detailed information on the geometry of the assembly and the basic operation of the reactor is provided. This 
information encompasses most of the inputs needed to simulate the assembly containing SF97. These inputs 
are summarized in Table [III However, there is insufficient information to simulate the full core, so by definition 
all Takahama simulations are assembly simulations. The fuel density used in simulations is an effective density 
to account for the variations in fuel pellet packing. This effective de nsity must be les s than 10.96 g/cm^, the 
theoretical density of UO2, but no more details are pro vided by Nakah ara et al.l 2001 1. For this reason, we use 
the value 10.07 g/cw? suggested bv lRoaue et al.l [2004 1. Another popular assumption is 95% of the theoretical 
density of UO2. 

The last important input to the simulation is the time dependent power of the reactor core. The type of 
power information provided by the Takahama benchmark is particularly unique. Simulations usually rely on 
the individual assembly power densities with the full core thermal power. The destructive analysis allows the 
use of the ^''^Nd concentration to produce a detailed power history along the rod, as displayed in Fig[T] The 
i48]\jj tech nique has a 3% uncertainty compared with uncertainties of <2% for the full core thermal power, 
iDiurcic et al.. (20Q9i |. The impact of this and other key uncertainties is discussed in the following section. 



5. Results for Rod SF97 

DRAGON uses the description of the fuel assembly from Table [llj to simulate the Takahama assembly. This 
simulation is performed in time steps with different power and duration as outlined in Fig. |TJ To simplify 
the model, we assume a constant temperature along the rod of 600 K instead of the detailed temperature 
gradient shown in Table HI We also assume a constant non-burnable concentration of boron in the moderator of 
630 ppm. To increase the speed of the calculation, two gadolinium rods are added to the assembly bringing the 
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Figure 2: The comparison of DRAGON simulation results to the measured mass inventories. Top: the DRAGON results 
are compared to other standard codes for reactor core simulation. Bottom: the critical inputs to the simulation are 
varied within their uncertainties. 



total to 16. This creates an eight-fold symmetry, and allows the simulation of a 1/8 segment of the assembly. 
A nuclear cross-sec tion library must be provided for the calculation. Our nominal simulation is performed 
using ENDF/B-VI, McLane et all [l996j . and we compare these results to those obtained with JENDL 3.2, 
iNakagawa et all |l995l |. 

The results of this DRAGON simulation are compared to the results from the destructive assay of the fuel 
rod. These are presented in Fig. [2] as a function of the sample's pos ition along the rod. The codes used for 
comparison range from ORIGEN 2.1, a simple fuel depletion model, iLudwig and Renieil Il989l. SCALE4 .4a, 
a code that uses an effective ID geometry to model the neutr on transport , DeHart and HermannI |l998{ . to 
the full Monte Carlo treat ment employed by MONTEBURNS, iDalle et al.l [2009l| . The code most similar to 
DRAGON is HELIOS. DcHart and HermannI f 1998j |. another 2D deterministic lattice code. 

We only present the results for ^^^U, ^'^^Pu, and ^^^Pu in Fig. [5] because the depletion ^ssy jg ^.j-^g game 
order as the uncertainty in the destructive assay, and therefore a useful comparison is not possible. For sample 
SF97-1, all codes show substantial deviations in the inventories of the plutonium isotopes. This sample is near 
the top of the rod. Therefore, the modeling of neutron leakage is difficult leading to these large deviations. 
Neglecting sample SF97-1, we calculate the average deviation along the rod. For ^'^^U, DRAGON's deviation 
is 3.2%. For ^"^^Pu, and ^^^Pu, the deviation is -1.3% and -3.9% respectively. For comparison the range for the 
other codes is -2.2% to 4.5% for ^ssy, -0.6% to 6.5% for ^aspu, and -4.0% to 3.4% for 24ipu. The DRAGON 
results have comparable accuracy to those from the other standard codes. 

Most work benchmarking reactor simulation codes stops here, but this neglects the uncertainty in these 
simulated results. To understand the effect of the input uncertainties on the simulation, we vary key input 
parameters within their uncertainties. We confirm that these key par ameters are the powe r, moderator temper- 
ature, boron concentration and fuel density, as originally outlined in lDiurcic et al.l |2009l |. We vary the power 
by 3%, the uncertainty in the -'^^^Nd method. We vary the moderator temperature by ±50 K, the variation 
along the length of the rod, and the boron concentration by 10%. The fuel density is varied by 1.5%, but we 
note that a deviation as large as 3% is needed to cover the range of values suggested for this effective density. 

The largest effect for ^ssy jg |-j^g power variation, a 6% effect. Since ^'^^Pu and ^^^Pu are the result of fast 
neutron reactions on ^^*U, these isotopes are sensitive to the total amount of fuel. Therefore, they are sensitive 
to the fuel density, a 2.5% effect. The results of these variations are summarized in the bottom of Figl^l If the 
larger 3% fuel density uncertainty is assumed, then the DRAGON results are consistent with both the measured 
mass inventories and the results of the other codes. We also find that the choice of cross-section library induces 
a 1.2% change in the ^35^ prediction, a 0.8% change for ^''^Pu and no change for ^'^^Pu. 

The Takahama benchmark focuses on the measured mass inventories. However, the quantities of interest for 
experiment are the fission rates, fi{t) in Eq.[Tl used to calculate the antineutrino flux. The mass inventories are 
proportional to the integrated number of fissions and the instantaneous fission rates. Because the systematic 
uncertainties on the Takahama inputs are large and the provided power is in a form not available for most 
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reactor cores, calculating systematic uncertainties on fission rates directly from the Takahama results is not 
possible. However, it is possible that an experiment could use these results to constrain the uncertainties in 
the flux calculation in its analysis. Regardless, this work has outlined a procedure for evaluating the major 
uncertainties and their effect on the fission rates. This work is applicable to all experiments that would use 
reactors as their source of antineutrinos. 



6. Conclusion 

The comparison of measured mass inventories from the Takahama-3 reactor to the DRAGON simulation 
demonstrates that the results obtained with DRAGON are accurate and of equal quality to other widely used 
codes for reactor modeling. We oiitline the procedure for evaluating the inpiit imccrtaintics and their effect on 
the fission rates that are critical to the antineutrino flux calculation. There are many interesting measurements 
to be done with reactor antineutrinos, and DRAGON is an excellent tool that should be used to aid these 
measurements. 
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